Mechanical Stability of a Strongly-Interacting Fermi Gas of Atoms 
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A strongly-attractive, two-component Fermi gas of atoms exhibits universal behavior and should 
be mechanically stable as a consequence of the quantum mechanical requirement of unitarity. This 
requirement limits the maximum attractive force to a value smaller than that of the outward Fermi 
pressure. To experimentally demonstrate this stability, we use all-optical methods to produce a 
highly degenerate, two-component gas of ®Li atoms in an applied magnetic field near a Feshbach 
resonance, where strong interactions are observed. We find that the gas is stable at densities far 
exceeding that predicted previously for the onset of mechanical instability. Further, we provide a 
temperature-corrected measurement of an important, universal, many-body parameter which deter- 
mines the stability — the mean field contribution to the chemical potential in units of the local Fermi 
energy. 



Strongly interacting Fermi systems are expected to ex- 
hibit universal behavior jlj]. In atomic gases, such strong 
forces can be produced in the vicinity of a Feshbach res- 
onance, where a bound molecular state in a closed exit 
channel is magnetically tuned into coincidence with the 
total energy of a pair of colliding particles 0|- In this 
case, the zero energy scattering length as, which charac- 
terizes the interactions at low temperature, can be tuned 
through ±oo. For very large values of \as\, the impor- 
tant properties of the system (e.g. the effective mean 
field potential, the collision rate, the superfluid transi- 
tion temperature, etc.) are predicted to lose their depen- 
dence on the magnitude and sign of as, and instead be- 
come proportional to the Fermi energy with different uni- 
versal proportionality constants. For this reason, table- 
top experiments with strongly interacting atomic Fermi 
gases can provide measurements that are relevant to all 
strongly interacting Fermi systems 0, 01 , thus impacting 
theories in intellectual disciplines outside atomic physics, 
including materials science and condensed matter physics 
(superconductivity), nuclear physics (nuclear matter), 
high-energy physics (effective theories of the strong in- 
teractions), and astrophysics (compact stellar objects). 

In a gas of degenerate fermions, compression of the gas 
is resisted by the Pauli exclusion principle, leading to an 
effective pressure known as the Fermi pressure or degen- 
eracy pressure. The Fermi pressure plays an important 
role in nature, providing, for example, the outward force 
which stabilizes neutron stars against gravitational col- 
lapse. Because of the Fermi pressure, a confined cloud 
of degenerate fermions is always larger than a cloud of 
bosons with an equivalent temperature and particle num- 
ber. This effect has been directly observed in two elegant 
experiments 



An important question is: At what point, if at all, do 
strong attractive interactions overcome the Fermi pres- 
sure? Previously, it was predicted that an atomic Fermi 
gas could become mechanically unstable for sufficiently 
large attractive interactions If true, this prediction 
presents a possible roadblock to attempts to use a Fes- 



hbach resonance to produce a high-temperature super- 
fiuid 0, S ' ttis should be unstable in the re- 
quired regime. Recently, however, it has become appar- 
ent that the gas may indeed be stable. Heiselberg has 
estimated the maximum ratio of the attractive poten- 
tial to the local Fermi energy for a strongly interacting 
Fermi gas. Using a self-consistent many-body approach, 
he finds that a two-component Fermi gas is mechanically 
stable, although a gas with more than two components is 
not 0. A similar conclusion can also be drawn from an 
examination of the compressibility of a strongly attrac- 
tive two-component gas [2| . 

In this Letter, we first show heuristically that the max- 
imum inward force arising from attractive interactions in 
a two-state gas of atomic fermions is limited by quantum 
mechanics to a value less than the outward Fermi pres- 
sure. We then demonstrate this idea experimentally by 
directly cooling a two-component Fermi gas of ^Li atoms 
to high degeneracy in an optical trap at magnetic fields 
near a Feshbach resonance, where strong interactions are 
observed ^^l- The highest densities obtained in the ex- 
periments far exceed those predicted previously 0] for 
the onset of mechanical instability. Imaging the cloud 
both in the trap and after abrupt release, we find no ev- 
idence of instability. To quantitatively describe the sta- 
bility, we present a measurement of the relevant universal 
many-body parameter (3 — the mean-field contribution to 
the chemical potential in units of the local Fermi energy. 
In Ref. |lflj |. we provided a first estimate of (3 assuming 
a zero temperature gas. Here, we present a method to 
properly account for the non-zero temperature of the gas. 
This revised measurement of /? is in good agreement with 
a prediction of 0|. 

The mechanical stability of the gas can be understood 
through a heuristic discussion of the forces acting on a gas 
of fermions in a 50-50 mixture of two spin components. 
The equation of state for a normal, zero-temperature 
Fermi gas is 



eF(x) + ?7mf(x) + ?7trap(x) ^ n, 



(1) 
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where /z is the chemical potential, eF(x) is the local 
Fermi energy, [/mf(x) is the mean field contribution 
to the chemical potential, and /7trap(x) is the trap po- 
tential. In the local density approximation, eF(x) — 
fc|,(x)/(2M), where kp is the local Fermi wavevec- 

tor, which is related to the density according to n(x) — 
fc|^(x)/(6^2) ^ 

For each spin component, there is an effective out- 
ward force arising from the Fermi pressure Ppcrmi = 
2n(x) eF(x)/5, even in the absence of interatomic interac- 
tions. The outward force per unit volume is — VPFcrmi- 
The corresponding force per particle is —(VPFcrmi)/"-- 
The local force per particle is then easily shown to be 



' Fcrr 



-VeF(x), 



(2) 



Since the density decreases with distance from the center, 
Fpermi is outward. 

Two body scattering interactions make a contribution 
Umf to the chemical potential, which is given in the mean 
field approximation by. 



Umf{x) = — — n(x). 



(3) 



where aos is an effective scattering length, which gen- 
erally depends on the local thermal average of a 
momentum-dependent scattering amplitude. The po- 
tential is attractive when Oog < 0, and repulsive when 
Oeff > 0. The local force per particle arising from the 
mean field contribution is just 



Fmf = -V/7mf(x). 



(4) 



If we assume that aoff is energy-independent and equal 
to as for as < 0, it is easy to show that the gas be- 
comes unstable for suitably large values of as- In this 
case, the inward force from the mean-field potential, 
Fmf oc Vn(x) while the outward force from the Fermi 
pressure, Fpermi oc Vn^/'^(x). Then the inward force 
exceeds the outward when |Fmf| > |FFermi|, i-e., when 
kp \aes\ > 7r/2. The corresponding density for mechani- 
cal instability satisfies n > no where no = 7r/(48|asp). 
This result matches the previous prediction of |^ which 
was derived via a more rigorous calculation. 

The assumption that Oeff is energy-independent, how- 
ever, is only valid if kp \as\ ^ 1. Outside that regime, 
this assumption violates the quantum mechanical re- 
quirement of unitarity. At intermediate densities 0, 



where as ^ kp^ 



cx n 



-1/3 



R, with R the range 



of the collision potential, two-body scattering is dom- 
inant, but the mean field interaction is proportional 
to a momentum-dependent two-body T-matrix element, 
which in turn is proportional to the scattering amplitude 
f{k). It is well known that / has a magnitude which is 
limited by unitarity to a maximum of 1/fc Hi. Hence, 
one expects that in a zero-temperature Fermi gas, one 
should use an effective scattering length in Eq. O with a 



maximum magnitude on the order of a^s = l/kp. For 
aoff = —1/kp, it is easy to show that the mean field con- 
tribution to the chemical potential is 



J7mf(x) = /3eF(x), 



(5) 



where the universal parameter /3 has a maximum value 
of — 4/(37r) = —0.42 in this heuristic treatment. In this 
limit, both the inward force, Fmf and the outward force, 
Fpcrmi SLTC proportional to the same power of the atomic 
density. The ratio of their magnitudes, however, is given 
by |FMF|/|FFcrmi| = < 1- Thus, one finds that a two- 
component Fermi gas is mechanically stable as a result of 
the quantum mechanical requirement of unitarity. This 
is in contrast to attractive Bose gases and Bose- 

Fermi mixtures J,T \ which exhibit dramatic instabilities. 

A possible criticism of our heuristic argument is that 
the two-body relative wavevector k lies in the range < 
k < kp. Since / increases as k decreases, f{k) should be 
averaged over Pr(k), the probability distribution for k, 
where / (ikVF(k) = 1. For a noninteracting gas at zero 
temperature, one can show 



W{k) = 
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Q{kp-k), (6) 



where k — |k|. Near a Feshbach resonance, where 
\as\ ^ R, we assume that f{k) takes the form expected 
for a zero energy resonance, and that the effective scat- 
tering length is determined by the real part of — /, i.e., 
aoff = (as/(l-l-fc^ where (...) denotes averaging with 
Eq.El We find that |aefr| has a maximum value of 1.05/ kp 
when kp \as\ = 2, and decreases slowly for A:i?|as| > 1, 
reaching 0.6/kp at kp\as\ = 10. 

Our heuristic estimate, which neglects the effects of the 
interactions on the free particle wavefunctions, yields a 
value of P comparable to that estimated by Heiselberg 0] . 
In contrast to our heuristic approach, the self consistent 
many-body approach of Ref. |1| also predicts that f3 is 
independent of the sign and magnitude of kp as in the 
intermediate density limit where /c^ jasj ^ 1. Hence, 
for strongly interacting fermions, /3 is always negative 
and the effective mean field interaction should always be 
attractive. 

To demonstrate that a strongly attractive, two spin- 
component Fermi gas of atoms is stable, we employ a SO- 
SO mixture of the two lowest hyperfine states of ^Li, i.e., 
the \F = 1/2, M = ±1/2) states in the low field basis. 
This mixture has a predicted broad Feshbach resonance 
at 860 G m A magnetic field of 910 G is apphed to 
produce a very large and negative zero-energy scattering 
length. 

The gas mixture is prepared and rapidly cooled to de- 
generacy at 910 G by forced evaporation in our ultra- 
stable CO o la ser trap, as described for our previous ex- 
periments [lO|. For our trap, lij±_ = 2n x (6625 ± 50 Hz), 
Lu^ = 2tt X (230 ± 20 Hz), and Cu = {Lu]_Lu^y/^ ^ 2tt x 
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(2160 ±65 Hz). Following evaporation, the trap is recom- 
pressed to full depth over 0.5 s, and allowed to remain 
for an additional 0.5 s, to ensure thermal equilibrium. 
The CO2 laser power is then extinguished and the gas 
is imaged as described in our previous paper ^f" 
ter release, the gas expands hydrodynamically, rapidly 
increasing in the transverse dimension while remaining 
nearly stationary in the axial direction 

We determine the number of atoms by numerically in- 
tegrating the column density (Table P). For a typical 
number, TV ~ 8.0 x 10^ atoms per state, the correspond- 
ing Fermi density for the experiments is calculated to 
be np — 4.8 x 10^^/cm^ per state 0. The density 
np, can be compared to the density uq = 7r/(48|a5p) 
predicted for the onset of mechanical instability, assum- 
ing a momentum-independent scattering length. Using 
the best available molecular potentials, which are con- 
strained by measurements of the zero crossing in the s- 
wave scattering length 0, li^l , the zero energy scatter- 
ing length as is estimated to be ~ — 10^*00 at 910 G 
(ao = 0.53 X 10^^ cm), within a factor of two 21]. Then 
no = 4.4x 10^^/cm'^, showing that the density np exceeds 
no by one to two orders of magnitude. 

Although the molecular potentials are constrained by 
the measured zero crossing in the scattering length 
lipj ]. the precise location of the Feshbach resonance still 
may be incorrect. Since the scattering length has not 
been directly measured, it is possible that the attractive 
potential is not as large as expected. However, a simple 
argument based on our previous experiments shows that 
as must be large ^3 ■ We observe similar hydrodynamic 
expansion after release from either full trap depth Uq or 
from a reduced trap depth of Uq/IOO. In the latter case, 
we estimate that kp^= 4000 oq. Since the observed hy- 
drodynamic expansion appears to be independent of the 
trap depth, as expected for a unitarity limited interac- 
tion, we conclude that in the shallow trap we must have 
jfeFisI > 1- Then we must have \as\ > 4000 ao and 
no < 6.9 X lO^Vcm^ 

Thus, we have clearly demonstrated that the gas is 
mechanically stable in the intermediate density regime, 
where \as\ is large compared to the interparticle spac- 
ing. The remainder of this paper provides a revised esti- 
mate of the universal parameter P which quantitatively 
determines the stability of the gas. In our previous esti- 
mate assumed a zero temperature gas and cal- 
culated (3 from an estimate of the release energy. There, 
we obtained a value of f3 substantially smaller in magni- 
tude than predicted 0. Here, we determine from the 
transverse spatial widths of the expanding gas and prop- 
erly include both finite temperature and hydrodynamic 
scaling effects. 

The one dimensional transverse spatial profiles after 
expansion for 0.4-0.8 ms are very well fit by normal- 
ized finite-temperature Thomas-Fermi (T-F) distribu- 
tions which determine the width cTx as well as the ra- 
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FIG. 1: Transverse spatial profile of the expanding cloud at 
600 ^s after release. 

tio of the temperature to the Fermi temperature, T/Tp. 
We expect that T/Tp is approximately constant during 
the expansion, since we observe hydrodynamic scaling of 
the transverse radii consistent with an effective poten- 
tial cx n^/'^ 01, suggesting an adiabatic process. The 
T-F shape is not unreasonable despite a potentially large 
mean field interaction, since the mean field contribution 
to the chemical potential, Eq. Elis proportional to the lo- 
cal Fermi energy. In this case, assuming Eg . [H for a nor- 
mal degenerate gas at zero temperature [llj , it is easy to 
show that the mean field should simply scale the Fermi 
energy of the trapped cloud without changing the shape 
from that of a T-F distribution 10]. The zero temper- 
ature spatial distribution then corresponds to that of a 
harmonic oscillator potential, with the frequencies scaled 
so that Lu[ — + P, i = x,y,z. Hydrodynamic ex- 

pansion then preserves the shape of the trapped atom 
spatial distribution 0, ^]. Assuming that T/Tp is 
small, one expects that a finite temperature spatial dis- 
tribution for a harmonic potential is a reasonable approx- 
imation. A Sommerfeld expansion [22l ] of the normalized 
density for \x\ < ax yields 

n{x)/N = -i^ [Ux] + 57r2 [T/Tpf h{x)\ , (7) 

where /o (x) = g^( x), f2{x) = g{x)/8 - g^{x)/6, and 
g{x) = y/l-x^/al. 

Eq. \7\ is used to fit the measured transverse spatial 
distributions. Fig. ^ We obtain the results given in Ta- 
ble^ for four trials each at expansion times t of 0.6 ms 
and 0.8 ms. Note that Eq.|7| begins to break down near 
\x\ ~ ax for T/Tp > 0.15. However, an exact treatment 
using polylogarithm functions yields similar results even 
for our highest temperatures where T/Tp ~ 0.18. 

The value of the parameter P can be determined 
from the transverse radii of the trapped cloud, ax{0) = 

^2e'p/AIuj'f, where e'p = hui' {6N)^^^ is the Fermi energy 
including the mean field contribution. Then, ax{0) = 
(1 -I- P)^/*axP 123, where axP = ^J2^7/{M1J\) , and 
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TABLE I: Value of /3: t is the expansion time, A*' is the number 
of atoms per state, a^p is Fermi radius without the mean field 
interaction, (t) is the measured Fermi radius obtained from 
the finite temperature fit to the data using Eq. |7] 





N 




T/Tf 


cr^it) (^m) 




600 


66,000 


3.49 


0.128 


99 


-0.237 




86,400 


3.65 


0.144 


102 


-0.220 




84,300 


3.63 


0.140 


101 


-0.234 




80,200 


3.60 


0.141 


101 


-0.208 


800 


67,200 


3.50 


0.146 


133 


-0.176 




87,000 


3.65 


0.179 


131 


-0.344 




70,400 


3.53 


0.151 


130 


-0.273 




82,000 


3.62 


0.183 


128 


-0.382 



ep = hCj {6Ny^^ is the Fermi energy in the absence 
of interactions. As we have pointed out previously, the 
observed anisotropic expansion can arise from unitarity- 
Umited colhsional hydrodynamics or superfluid hydro- 
dynamics In either case, the effective potential is 
cx n^/'^, and we can assume ax{t) — <Jx{0)bx{t) [lollll|. 
Hence, 

/5=fr^)'-l, (8) 

For our trap parameters, bx{0.6ms) = 29.74 and 
6a; (0.8ms) — 39.88. The scale factor bx{t) properly in- 
cludes the spatial anisotropy of the expansion and the 
correct hydrodynamic scaling. We obtain from Table Q 
an average value (3 = —0.26 ± 0.07. 

This result is consistent with that obtained from the 
measured transverse release energy. The release energy 
can be initially calculated from the zero temperature T-F 
fits |lOj or calculated numerically from second moment 
of the density. Since T/Tp ~ 0.15 from TableHl the mea- 
sured release energy is larger than the zero temperature 
release energy by a factor ~ l-h(27rV3){T/TF)2 = 1.15. 
Reducing the release energy by a factor rj and using the 
method of Ref. ^3 > yields an average value of consis- 
tent with that given above. 

The measured average value of = —0.26 ± 0.07 can 
be compared with predictions of the energy per particle 
for the conditions of our experiment, where kpas — —7.4. 
For our equation of state, the sum of the local kinetic and 
mean field energy per particle is (3/5)(l -I- /3)eF(x) [23| . 
In Ref. 0, Heiselberg provides two estimates for the to- 
tal energy per particle. His Eq. 1J_, based on the Galitski 
equations, yields /?caic = —0.54 |23. Using the Wigner- 
Seitz cell approximation, his Eq. 14 provides an alterna- 
tive estimate, /Scale — —0.33. An independent calculation 
by Steele using effective field theory yields (3 — —0.46. 
The theoretical predictions are for zero temperature, and 
the second is in reasonable agreement with our measure- 
ments including the temperature correction to the spatial 



distributions. However, the predictions do not include 
the additional temperature dependence of (3 which arises 
from a thermal average. Since this important universal 
many-body parameter can now be experimentally mea- 
sured, further refinement of both the experimental mea- 
surements and the theory is worthwhile and may permit 
the observation of the superfluid corrections to (3. In 
addition, mixtures of the three lowest hyperfine states 
of ^Li may permit observation of predicted mechanical 
instabilities in a three-component Fermi gas 0|. 
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